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Abstract 

We formulate the three dimensional Thirring model on a spacetime lattice and 
study it for various even numbers of fermion flavors Nf hj Monte Carlo simulation. 
We And clear evidence for spontaneous chiral symmetry breaking at strong coupling, 
contradicting the predictions of the 1/Nf expansion. The critical point appears to 
correspond to an ultra-violet flxed point of the renormalisation group; a fit to a RG- 
inspired equation of state in the vicinity of the flxed point yields distinct critical 
exponents for Nf = 2 and Nf = 4, while no fit is found for Nf = 6, suggesting 
there is a critical number Nfc < Q beyond which no chiral symmetry breaking occurs. 
The spectrum of the Nf = 2 theory is studied; the states examined vary sharply but 
continuously across the transition. 
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1 Introduction 


The Thirring model is a quantum field theory of fermions interacting via a current- 
current contact term, described in three dimensional continuum Euclidean spacetime 
by the following Lagrangian: 

L = + m)i)i + ( 1 . 1 ) 

Here 'ip and ip are four-component spinors, m is a bare mass, and i runs over Nf 
fermion species. After the introduction of an auxiliary vector field, the above La¬ 
grangian can be rewritten as: 

TV. ^ 

^aux = +i4 + (1-2) 

which now contains the fermionic bilinear term of the usual QED lagrangian but 
with a mass term for the vector field, which spoils the gauge-invariance of Caux- 
This is no surprise, since the original lagrangian did not show any local symmetry. 
The global symmetries are a U{Nf) one, corresponding to rotations in flavor space, 
and, for m —>■ 0, a U{Nf)A_ chiral symmetry. Since the coupling constant g has mass 
dimension [g] = (1 —|) = —the usual perturbative expansion is non-renormalisable. 
However perturbation theory is not the end of the story and, as suggested several years 
ago [|, 1^, different approaches can yield a well-defined continuum limit, corresponding 
to a UV-stable fixed point of the renormalisation group (RG). 

The Gross-Neveu model (GN), which has been studied both numerically and in 
the 1/A^/ expansion 0, Q, exhibits this behaviour. The model turns out to be 
renormalisable to all orders in 1/A^/, with a UV-stable RG fixed point at guv 4 0 
and chiral symmetry breaking for g > guv] only for g ~ guv can the ratio of the 
physical scale (either a particle mass or a scattering cross-section) to the UV cutoff 
be made small. The theory defined at this point describes the IR dynamics of a linear 
sigma model with a Yukawa coupling to fermions Q, which is super-renormalisable 
about the origin: 

I^Yu = V’(^ + w -t- m)pj + ^(^0-)^ + (1-3) 

For a critical value of v in the IR regime the kinetic and quartic terms of the a-field 
become irrelevant operators and thus the system equivalent to the GN model. The 
critical behaviour for the GN model can be summarised by the behaviour of the [3- 
function, as shown in Fig. ^ In two dimensions, the theory is asymptotically free. 
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and the /^-function exhibits the usual UV hxed point at the origin. For 2 < d < 4, 
the fixed point is shifted away from the origin, separating a strong coupling regime, 
where dynamical chiral symmetry breaking is possible, from a weak coupling phase 
where the theory is chirally symmetric |^. 

The situation for the Thirring model is less clear. At both leading and next-to- 
leading order in 1/Nf the renormalisability of the massless model has been estab¬ 
lished p. For d < 4 the vacuum polarisation is actually UV-finite so long as the 
regularisation respects current conservation. As a consequence a continuum limit 
may be defined for any value of the coupling g, which would control, e.g., the ratio 
of physical fermion mass to vector bound state mass. In other words, the /3-function 
is vanishing for any value of g and the interaction is a marginal operator, 

whereas the interaction for the GN model is relevant 0. Another prediction of the 
1/Nf expansion is the vanishing of the vector bound state mass in the strong cou¬ 
pling limit mg —> cx). A massless particle mediating an interaction between conserved 
currents suggests an equivalence between the UV behaviour of Thirring model and 
the IR one of QEDg, which would thus be the analogue of the equivalence between 
the GN and Yukawa models. Such a conjecture is supported by the fact that the 
0{1/Nf) corrections in the two models appear to coincide [^, |^. Moreover in the 
1/Nf expansion, as in perturbative gauge theory, no diagram ever appears resulting in 
the dynamical breaking of chiral symmetry leading to spontaneous mass generation. 

The prediction of no spontaneous mass generation is contradicted, however, by 
results from a completely non-perturbative approach, namely solution of Schwinger- 
Dyson (SD) equations, which predicts a non-vanishing chiral condensate {‘ip'ip) in the 
limit m —>■ 0 0.0.0. In this approach a non-trivial solution for the dressed fermion 
propagator S{p) = -|-S(p))“^ is sought, i.e. one in which the self-energy S(p) 

and hence ('0'0) are non-vanishing in the chiral limit m —> 0. Unfortunately, the SD 
equations can only be solved by truncating them in a somewhat arbitrary fashion. 
The usual approximation p, |T^, is to assume that the vector propagator is given by 
its leading-order form for m = 0 in the 1/Y/ expansion, viz. 





k^,ku 

k‘^)\ 8 j k^ 


(1.4) 


and that the fermion-vector vertex function is well-approximated by the bare vertex 
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(the so-called “planar” or “ladder” approximation); 




w 


=7m- 


(1.5) 


The most systematic treatment has been given by Itoh et al [0, who note the 
eqnivalence of (|1.2|) to a gange-hxed form of a ganged fermion-scalar model. This 


allows the specification of an alternative non-local gange-hxing condition, which in 
tnrn entails A{p) = 1, with no effect, of conrse, on gange-invariant qnantities snch as 
the physical mass or {'ip'ip)- The resnlting SD eqnations can be solved exactly in the 
limit > oo; a non-trivial solntion for S(p) exists for 


N) < ~ 4.32. 


( 1 . 6 ) 


This is identical to the critical Nfc predicted for non-trivial IR behavionr in QEDg 
r3| . However, since the integral eqnations reqnire the introdnction of a UV cntoff A, 


a featnre of this solntion is that the indnced physical scale fi depends on Nf in an 
essentially singnlar way: 


/ 


p 
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cx exp 


27r 
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A££ _ 1 

Nf 


_ 13 

(i/ji/j) cx A 2 /i 2 cc exp 


Svr 




V 


Nfc _ 4 
Nf 


(1.7) 


This implies that in the strong conpling limit a continnnm limit only exists as Nf ^ 
Nfc 

The behavionr o is an example of an infinite order, or “conformal” phase 
transition, originally discnssed in a particle physics context by Miranskii and collab¬ 
orators in qnenched QED 4 [n|, and recently discnssed in more generality HI- One 
prediction of [^] is that althongh the order parameter (' 0 ' 0 ) is continnons at the 
transition, there is a discontinnity in the spectrnm of light excitations, with no light 
scalar resonances on the symmetric side of the transition. In the broken phase, of 
conrse, Goldstone’s theorem predicts light pions. This behavionr has recently been 
exhibited in an investigation of fonr dimensional SU(A'"c) gange theory with an inter¬ 


mediate nnmber Nf of fermion flavors a solntion of the Bethe-Salpeter eqnation 
in ladder approximation fonnd no low-momentnm pole in the symmetric phase. 

Unfortnnately no analytic solntion to the SD eqnation for the Thirring model 
exists for g‘^ < 00 ; however nsing different techniqnes Kondo |]^ has argned that 
there is a critical line Nfc{g‘^) in the {g"^, Nf) plane, which is a smooth invertible 
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function. Therefore for integer Nf < Nf^ one might conjecture the form 

( \ 

_ n 

('0V’) oc exp 


V 


31-1 

n2 


( 1 . 8 ) 


corresponding to a symmetry restoring transition at some critical point = g‘^{Nf). 
Presumably in this scenario the coupling g has become relevant: there may exist 
a novel strongly-coupled continuum limit at the critical point not described by the 
1/A^/ expansion, which is simultaneously the UV limit of the Thirring model and the 
IR limit of QEDg. Since there is no small dimensionless parameter in play, this hxed 
point would be inherently non-perturbative. It is worth mentioning in passing that 
a novel IR hxed point in QEDg has been proposed by Aitchison and Mavromatos in 
connection with non-Fermi liquid behaviour in the normal phase of superconducting 
cuprates |T^ . 


Although this scenario of a non-trivial hxed point for Nf < Nfc is attractive, there 
are good reasons to be cautious. Using a diherent sequence of truncations Hong and 
Park have found chiral symmetry breaking for all Nf |j^, with 


1 

— oc exp 
9c 


Nf^ 

16 


( 1 . 9 ) 


a result which is also non-analytic in 1/A^/ and hence beyond the reach of the 1/A^/ 
perturbative series. Moreover in QEDg, studies beyond the planar approximation, 
using improved ansatze for the vertex P^, suggest that the condition A{p) = 1 is 
unphysical, and that chiral symmetry is spontaneously broken for all Nf In the 

current context this would imply g1 < oo for all Nf. 

Because of the systematic uncertainty attached to the SD approach, and because of 
the intrinsic interest in exploring the non-perturbative behaviour of a fermionic model, 
we have studied the three dimensional Thirring model by Monte Carlo simulation on a 
spacetime lattice. Our initial results appeared in |]^, where we claimed that critical 
points with distinct critical indices exist for the models with Nf = 2,4, but that 
evidence for a critical theory with Nf = 6 was not found. In this paper we summarise 
those results, and present further data which enable us to rehne our description of 
the critical behaviour of the Nf = 2 model. In Sec. 2 we give a fairly lengthy 
exposition of the lattice formulation, showing that the lattice method is itself not 
without systematic uncertainties. Although we use staggered lattice fermions which 
have a continuous remnant of the U(A^/)a axial symmetry, a transcription to a basis 
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in which the Fermi helds have explicit spin and flavor indices reveals that there are 
many other interactions present at tree level besides the desired current-cnrrent form. 
Perhaps a robust counter-argument to this problem is that for a strongly-coupled 
fixed point a non-perturbative regularisation is always necessary, and therefore “what 
you see is what you get”, i.e. the fixed point theory defines itself via its critical 
exponents and spectral quantities. In addition, the lattice form of the current is not 
exactly conserved, resulting in an additive renormalisation of the inverse coupling 
1/g'^. This makes it difficult in practice to identify the strong coupling limit —*• oo. 
A description of the simulation, and a list of the measured observables, namely the 
condensate, susceptibilities and propagators, follows. 

In Sec. 3 we present results for the equation of state, i.e. m = f {{'ijj'ip), g). Since it 
is impossible to work directly in the chiral limit m = 0 , we find it convenient to fit to 
a phenomenological form inspired by the standard power-law scaling form used in the 
theory of ferromagnetism, with its attendant critical exponents /?, 7 , 6 and u, together 
with the additional constraints of hyperscaling, and, from the ladder approach to the 
gauged Nambu - Jona-Lasinio model that 6 — 1/(3 = 'j = 1. We justify this 
approach not so much by any theoretical argument but because its shortcomings if 
any should manifest themselves as variations in the fitted values of the exponents as 
data is taken from systems on larger volumes and closer to the chiral limit. In ||^ 


the fitted exponents for the Nf = 2,4 models were distinct, a feature not seen in 
similar studies of QED4 


Here we supplement the data of []^ for Nf = 2 
with results from a larger lattice (16^) and a smaller bare mass m, and refine the fit 
by combining data from several lattice volumes in a finite size scaling analysis. We 
find critical exponents and g^ completely consistent with those of our earlier studies, 
which strengthens the case for a critical point described by power-law scaling. We 
also have some new results for the Nf = 6 model on 12 ^ systems, supporting our 
previous claim that there is no transition in this case, and hence that Nfc < Q for the 
Thirring model. 

In Sec. 4 we present exploratory results for susceptibilities, and for fermion and 
bound-state masses. Our lattice size of 16^ is unfortunately not long enough in the 
temporal direction to permit real accuracy; nonetheless qualitative features of the level 
ordering in either phase do emerge. The susceptibilities reflect the phase transition 
very nicely, with longitudinal and transverse susceptibilities xi Xt varying from 
near-degeneracy in the symmetric phase (thus indirectly supporting the hypothesis 
5 — 1//? = 1) to Xi S> Xz in the broken phase. The effects of dynamical fermion loops 
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in the quantum vacuum can also be estimated, and appear to become dominant 
in the broken phase. The fermion, pseudoscalar and scalar masses all behave as 
expected, the pion mass staying small and roughly constant, while fermion and scalar 
masses change sharply but continuously across the transition. There is no sign of the 
discontinuities described in [|^. Finally studies in the vector-like channel reveal light 
states both with vector and axial vector quantum numbers in the symmetric phase. 
This is interesting because no light axial vector is predicted in the continuum 1/A^/ 
expansion at leading order, as we show in Sec. 4. Both states become difficult to 
observe in the broken phase. Our conclusions are presented in Sec. 5. 


2 Lattice formulation 


We want to discretise the action of Eq. ( |1.1|) on a hypercubic Euclidean lattice. In 
order to recover a partial U{N)a axial symmetry in the massless limit, we use N 
flavors of staggered fermions. The lattice action then reads: 

^ = \Y. M^)Vi^{x)[Xi{x +fi)-Xi{x-fi)] + 

xfii 

mJ2xi{x)Xi{x) + Xi{x)Xi{x + fi)Xj{x + fi)Xj{x) (2.1) 

xi xfiij 


where y, y are the staggered fermion helds, the Kawamoto-Smit phases, m is the 
bare fermion mass, the flavor indices i^j run from 1 to and the vector character of 
the interaction is encoded in the geometrical structure of the four fermion term. 

To express the action as a bilinear in y, y, which makes it amenable to simulation, 
we introduce a bosonic auxiliary held as in (|1.2|) . In this paper we choose a form 
which we will refer to as non-compact 


^ + fi)+h.c. 

xfii 


N 

+m ^ Xi{x)Xi{x) + ^ 41 ^ 

xi ^9 xfj. 




N 


xyi J 


Xfl 


( 2 . 2 ) 


where we have introduced M[A,m] to denote the fermionic bilinear, which depends 
on both the auxiliary held and the bare mass. It is straightforward to integrate over 
the auxiliary held and recover Eq. 0. 
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In the chiral limit m —> 0 both the original action (^) and the bosonised form 
( p.2D are invariant under a global axial rotation, which for iV = 1 is written 


M exp(ia8{x))x(x) ; x(x) i-^ exp(!a£(i))\'(i), 


(2.3) 


where the phase e{x) = i)*i+*2+a;3_ generalisation to U(N) is straightforward. 

We will refer to this as the chiral symmetry of the model. In the m —>■ 0 limit for 
certain values of the coupling g‘^ > g‘l{N) it is spontaneously broken, signalled by the 
appearance of a condensate {xx) 7^ 0- 


Some comments on the lattice formulation are appropriate. Eq. (^) in the long 
wavelength limit (i.e. assuming ad'ip -C 'ip) is clearly closely related to the continuum 
auxiliary form (|1.2|) , and as such is suitable for a perturbative expansion in small-mode 
fluctuations of the helds, corresponding to the continuum 1/A^/ expansion. This 
expansion is well-behaved because to leading order the vector propagator = 

( 24 ^( 0 )24i,(x)) receives a contribution not just from the quadratic term but also 
from the vacuum polarisation loop 11^^,, which is hnite in d < 4. We hnd in the deep 
Euclidean region [|] (Cf. (|L^) ) 

1 




lim Dni,{k) (X V^v{k)^= + longitudinal 




where Vn^ is the transverse projector 




/if/ 


(*;) = 




(2.4) 


(2.5) 


Because the interaction current tpXij.'ip is conserved, the longitudinal component of 
has no physical signihcance [|I||T3|^. The 1/k behaviour of in effect suppresses 
short wavelength excitations of the auxiliary held, and justihes the long-wavelength 
expansion a posteriori. It means we might hope the form (|2.2| ) has the correct con¬ 
tinuum limit. This assumption may not be justihed, however, due to the chiral phase 
transition at hnite g^. Since the I/IV/ expansion predicts no phase transition, we 
know that there will be at least some part of coupling constant space where the 


above arguments cannot hold. In this case it is not quite so clear that (|2.2|) coincides 
with the continuum form ('ll . 1|) . 

To see why there might be a problem, we must consider the continuum havor 
interpretation of the lattice action (p.l|) . First dehne a unitary transformation to 
helds q, q 

1 


qriy) = 


Ui 

dr 


{y) 

{y) 


4^2? 


^ P 


Br 


Xi{P-,y), 


( 2 . 6 ) 
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where y denotes a site on a lattice of spacing 2a and P is a 3-vector with entries 
either 0 or 1 ranging over corners of the elementary cube associated with ?/, so that 
each side of the original lattice corresponds to a unique choice of {P]y). The 2x2 
matrices Tp and Bp are defined by 

Tp = T['Ti‘Ti> ; = (2.7) 


where are the Pauli matrices. It is then possible to recast the bilinear term of ( |2.1| ) 
as follows; 


^kin (2®) / A Qiiy) [Iu ® II 2 ) . T 

yfj. V 

m{2af'^qi{y){li®l2)qi{y). ( 2 . 8 ) 

y 

The direct product is between 4x4 Dirac matrices, acting on spinor degrees of 
freedom, defined by 


7/7 = 






(2.9) 


and 2x2 matrices acting on flavor degrees of freedom. Each lattice fermion flavor in 
three dimensions corresponds to two continuum four-component spinor flavors. Hence 


Nf = 2N. (2.10) 

In terms of this new basis the chiral symmetry of the model (|2.3|) reads 

g exp(iQ((75 0 l))g ; g gexp(ia;(75 0 1)). (2.11) 

The condensate {xx) fhe new basis reads (g(l 0 l)g). 

The hrst and third terms of (|2.8| ) are clearly of the same form as the continuum 
action for free massive fermions, whereas the second term, which is flavor non-singlet 
and Lorentz non-covariant, is formally 0{a) and hence hopefully irrelevant. The 
problem comes when we consider the interaction term. In the auxiliary held formal¬ 
ism, only in the long-wavelength limit would we expect the interaction to assume its 
expected continuum form A^qiipfyL ® ^ 2 )qi HH- When A is integrated out to recover 
the original form (p^.l|) though, all momentum modes are summed over with equal 
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weight. In fact, in four-fermi form the interaction contains many extra termsQ; 
^2 (2a)3 

- A 

ym 

n \ / \ _/\/. _ j-Xr- /\\^ 




9 J 

2N 2 ^ 

{9i{y){'ly ® 1)(1 + //.)?*(!/) + 9i{y){i'lA ® 

((1 + fyi)9ji.y)i.ly ® ^)9j{.y) + ^^l9j{y){.-^lA ® ^1)9: 

+ (%(|/)(757m < 8 ) 1)(1 + fMi(y) + %(2/)(*7574 ® 

((1 + fy)9jiy)ilyl^ ® ^)9j{y) + ^M?i(l/)(-*7475 ®r*^)<l: 

+ Y.(9i{.y){il4.1blul^, ® 0(1 + f^,)qi{,y) + qi{y){,lyl^ ® T*Tp6f,q, 

V 

((1 + /M)?i(2/)(0757M7!. ® O^jO) + ^iAAy){-iui^ ® ^*^^u)q3 

+ Z!(^i(2/)(“07i^7M ® ^u){^ + fi^)qi{y) + qi{y){-ii^ ® r*Tp6^qi 

V 

((1 + /M)?i(2/)(07/.7^ ® -rDqAy) + ^iAAy)^-^y ® T'mO^: 
where the difference operators / and 6 are dehned by 


X 


X 


X 


( 2 . 12 ) 


fMv) +f‘) + 'Av-IJ-)) ; S^q{y)^^{<i{y + ii.)-q(y-j£]). (2.13) 

The terms containing 5 are formally 0(a) and hence probably irrelevant, but no 
such argument can be used for the remainder, which, in addition to the desired form 
(%( 7 ai <81 l)q'j)^, contains extra flavor non-singlet and non-covariant terms. If the 1/iV/ 
expansion were valid, then one could argue that only interactions between conserved 
currents would generate long-range correlations between the A helds, and that all 
other interactions would become irrelevant in the continuum limit - a similar effect 
appears to be the case in other two [pTj ^ and three 0 dimensional four-fermi models 
where a renormalisable expansion is available. Even in this case there might also be 
contributions in the chiral limit from the interaction (^ 1 ( 7^75 0 However, as 

argued above, in the present case we know the 1/iV/ expansion fails qualitatively in 
the chirally broken phase. It therefore becomes a matter of experiment, by studying 
the current-current and AA correlations in various channels numerically, to determine 


the continuum limit of the action (| 2 . 1 |) . 

A second issue we must address is that despite the above comments there is an 
important difference between the continuum 1/A^/ expansion and the 1/N expansion 
^The corresponding result in four dimensions was first communicated to us by M. Gockeler pffl 
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applied to the non-compact lattice action (|2.2|) . This concerns vector cnrrent conser¬ 
vation. In the continnnm model a regnlarisation snch as Panli-Villars, which conserves 
the interaction cnrrent, must be assumed in order to ensure that the superhcial 0(A) 
contribution to vacuum polarisation vanishes, leaving a hnite result BBi- Now, in 
non-compact lattice QEDg, the vacuum polarisation is again UV hnite, but this time 
because of a cancellation between 0{a~^) contributions from two separate diagrams 
(Fig. see eg. [^). The second diagram appears because the gauge-invariant inter¬ 
action term in QED is xexp(ieA^)x, and hence contains n-photon-electron-positron 
vertices. In the action (|2.2D , on the other hand, the interaction is x'^A^Xj so only the 
left-hand diagram is present. Therefore cancellation of 0{a~^) divergences does not 
occur. We hnd for the vacuum polarisation tensor to leading order in 1/A^: 



{k) = (2.14) 

CL 


CONT 


where 



( 2 . 15 ) 


and, recalling a factor of two from ( p.lO|) 



(2.16) 


This means that for the inverse auxiliary propagator (setting a = 1) 


= ( 2 . 17 ) 


LATT 


CONT 


i.e. 



( 2 . 18 ) 


If we compare this with the continuum result 



( 2 . 19 ) 


we see that the 0{a ^) divergence can be absorbed by a wavefunction renormalisation 
together with a coupling constant renormalisation: 
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Now, the physics described by continuum 1/iVj perturbation theory occurs for the 
range of couplings G [0, cx)), i.e. for G [ 0 , 5 '^*^), where to leading order 


lim 

m —>0 


9ii 


2 

3' 


( 2 . 21 ) 


We therefore expect discontinuous behaviour at bare coupling g^ ~ gf^^, since at this 
point the auxiliary propagator is not defined (i.e. D~^ in ( p.lTp has no inverse). For 
9^ > 9fim^ since U^ont q in the limit —>• 0, becomes negative, suggesting 
that the model is not unitary in this region. 

In Fig. 1^ we plot a = {xx) vs. l/g'^ for N = 1,2,3 and m = 0.10 on a 8^ 
lattice 0. The three models have apparently coincident condensates in the strong 
coupling region l/g'^ < 0.3, which appear to be tending to zero in the g'^ ^ oo limit. 
It is tempting to associate this behaviour with g"^ > gf^^, but the correspondence with 
the value (|2.21|) is not good. It may well be that the numerical value of the right-hand 
diagram of Fig. is considerably reduced in a chirally broken vacuum; in any case we 
do not expect the 1/N expansion to be accurate here for the reasons described above. 

The ill-defined behaviour of the g‘^ ^ oo limit of the non-compact action (|2.2| ) 
presents no problem for most of the simulation results presented in this paper, which 
focus on the critical region of the N = 1 model around l/^f^ ~ 2.0, and of the N = 2 
model with l/gf^ ~ 0.6. However, we expect g^ to increase as N increases flilBio, 
ISl , so that as models with larger N are explored, eventually the relation between g"^ 
and gj^ will become important, particularly if lattice studies are to resolve the issue 
of the existence and value of Nfc- It is therefore worth making remarks on two other 
possible forms for the lattice Thirring model. 

For N = 1 there is an alternative but equivalent way of introducing the auxiliary 
held, which we will call the compact form: 


S = \Y.Xiix)vAx) U + 

xfii \ 


\2_l_ 

N 


exp{i9^{x)) \ Xi{x + fi) + h.c. 


( 2 . 22 ) 


where the link variables 0fj,{x) are freely integrated over |^, It is straightforward 
to perform this integration to yield a term cx —2g^x{x)x{x + fi)x{x + fi)x{x)), 
with If) a modihed Bessel function. If we dehne Jq by its power series expansion, then 
the Grassmann nature of the y, x truncates the expansion and forces equivalence 
with the partition function based on (O). We have performed some trial simulations 
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using (2 .22 ), and found identical results to those of the non-compact form (|2.2|) within 
statistical accuracy. It is interesting to note in passing that using a similar method 
one can show the equivalence of (|2.1| ) with the strong gauge-coupling (i.e. /3 = 0) 
limit of the gauge-fermion-scalar (xUcp) model recently studied in connection with 
non-trivial continuum limits in 3 and 4 dimensions 0. 


The bilinear action ( ^.221) can be extended to iV > 1 simply by introducing explicit 
flavor indices i on y, y. Now, however, terms of the form (x(a;)x(a; -|- /i)x(a; -|- fi)x{x))^ 
with n > 1 are not suppressed by the Grassmann nature of y: for N flavors of 
staggered fermion there are thus n-point couplings in the action with n < AN. Since 
ultimately we wish to understand and contrast the behaviour of the Thirring model 
for various values of we prefer the non-compact form (|2.2|) , where the interaction 
is of the same form for all N. It may well prove to be the case that in the RG sense 
the higher-point couplings are irrelevant, and hence ( p.2|) and (p.22|) have identical 
continuum limits. Both variants deserve study. 

Finally we mention a recent numerical study of a gauge-invariant version of the 
three dimensional Thirring model by Kim and Kim ||^; the action is 

^ = \Y.N{x)7]^{x)e^Y>{id^{x))xi{x +fi)+\i.c. 


X^l 


N 


+m Xiix)xi{x) - cos(0(x + /i) - 0(a;) + 0f,{x)), (2.23) 

xi 9 xfj. 

with (j) a gauge-covariant scalar held. This action has a local gauge invariance by 
construction, corresponding to the hidden local symmetry noted in Refs. ||I^ , p!7| . In 
unitary gauge 0=1 this action resembles (p.2|) in the long-wavelength limit, but now 
has a periodic dependence on the auxiliary held. The main result of is that the 
Nf = 6 model appears qualitatively diherent from Nj = 2, once again supporting a 
value of Nfc < 6. It is an open question whether the action ( 2.23 ) lies in the same 
universality class as (ED- 

We simulated dynamical fermions using a hybrid Monte Garlo (HMG) algorithm. 
In all cases we used a random trajectory length drawn from a Poisson distribution 
with mean 0.9, and adjusted the timestep to maintain an acceptance rate of 70% 
or greater. Typically we performed measurements on conhgurations separated by 
between 2 and 5 trajectories. Our estimated autocorrelation times measured on se¬ 
quences of conhgurations varied from ~ 1 to ~ 3 near the critical coupling; all our 
quoted errors take autocorrelations into account. The convergence criterion we used 
in the conjugate gradient matrix inversion was that the modulus of the residual vec- 


13 

























tor was 10“® per lattice site during guidance and 10“® per site for Metropolis and 
measurement. Unlike the case of the GN model P], we were unable to boost accep¬ 
tance rates by tuning the parameters of the guidance Hamiltonian. However, the 
vectorlike nature of the interaction allows the simulation of the model for any value 
of the number of lattice species N, corresponding to a number of physical fermions 
Nf = 2N, via an even-odd partitioning. For the very same reason, the only diagonal 
term in the fermion matrix comes from the mass term in the action. As m gets small, 
we hnd that the fermion matrix becomes very difficult to invert. Simulations with 
small m are therefore extremely time-consuming and represent a major problem for 
the extrapolation of the data to the chiral limit. Most of the results reported in this 
paper use between 200 and 500 measurements at each value of the coupling and bare 
mass. We have used periodic boundary conditions for the auxiliary held, and for 
the fermions used periodic in the spacelike directions and antiperiodic in the timelike 
direction. Most of the new results in the paper are obtained on 16^ lattices; however 
we also did some runs on an asymmetric 8^ x 16 lattice, to give some estimate of hnite 
volume effects on correlation function measurements. 

In our simulations we measured the following observables: 

(i) The chiral condensate 

= ixx) = (2.24) 

where M is the fermion kinetic operator introduced in eq.( |2.2| ). In Sec. 3 we will 
discuss numerical hts to the model’s equation of state using these measurements. To 
enable alternative htting procedures we list all our measurements made on symmet¬ 
rical (i.e. L^) lattices for Nj = 2 in Tabs. and ^ for Nf = A in Tab. ^ and for 
= 6 in Tab. 

(%%) The longitudinal susceptibility 

Xi = hn^xxi^) xx{x))c 

^ X 

= ^ ^[((trM-i)( trM-^)) - (trM-i)^] - (2.25) 

^ X 

— Xls T Xlnsi 

where we have denoted the components with disconnected and connected fermion 
lines respectively as “singlet” xis and “non-singlet” xins in analogy with QCD havor 
dynamics. 


14 


(in) The transverse susceptibility 


Xt 


1 

77lZ(X£X(0) X£X{x))c 

^ X 

Xts T Xtns 



(2.26) 


where we have used the fact that {x^x) vanishes, and the last line follows from a 
Ward identity. 

(iv) The “vector susceptibility” 


Xv 


( T. (n(o)nw), 

^ X,fM=l,2 
Xvs T Xvns: 


(2.27) 


where is the interaction current: 

= ■^{x{x)Vt^{x)x{x + h) + x{x + fl)v^,{x)x{x))■ (2.28) 

Due to the inhomogeneous boundary conditions, we only averaged this quantity over 
the spacelike directions. 

All susceptibilities, of course, are identical to the bound state propagator in the 
appropriate channel evaluated at zero momentum. The non-singlet components were 
therefore easily evaluated from the bound state correlators discussed below. The 
singlet parts, involving traces over a single power of M~^, were evaluated using gaus- 
sian noise vectors. Since the susceptibility signals are intrinsically noisy, we used 5 
noise vectors per conhguration, thus obtaining 5 independent estimates of (yy), and 
5 X 4/2 = 10 independent susceptibility estimates. 

In addition to these operators we have also performed spectroscopy by evaluating 
timeslice propagators. We examined the following channels: 

(v) Fermion propagator 

Cf{t) = Re ixi^) X{x,y,t)), (2.29) 

x^y even 

where note that the signal is improved if only sites an even number of translations in 
the transverse directions are included in the sum. 

(vi) Local composite operators 

Ci{t) = II(xrx(o) x'^x{x,y,t)), (2.30) 

x,y 
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where the space-dependent phase T is chosen to project onto various channels as 
follows: 


“pion” 

r(a;) = e{x), 

(2.31) 

“scalar” 

r(a;) = 1 , 

(2.32) 

local vector” 

r(a:) = (-1)"+* + (-l)^+h 

(2.33) 


The local vector operator is chosen by analogy with the rho-meson operator commonly 
used in QCD simulations. 

(vii) One-link composite operators 

“conserved vector” Ccv{t) = ^ XI H (2-34) 

^ At=i>2 

i.e. the two-point correlator between interaction currents Vji(x), corresponding to the 
non-singlet component of the vector susceptibility 

Of course, before assigning names to the various channels we should examine the 
spin/flavor structure of the fermion bilinears using the transformation (p.6|) . Since 
the staggered fermion action is invariant only under translations by an even number 
of lattice spacings, in general, as t ^ oo, two-point correlators behave as 

C{t) = (2.35) 

i.e. displaying signals in both “direct” and “alternating” channels, characterised by 
distinct masses Md and Ma respectively. The spin/flavor assignments we find are 
shown in Tab. It is interesting to note that in three Euclidean dimensions the 
local rho meson operator, which we have referred to as “local vector”, does not in 
fact have vector quantum numbers. 

In the calculation of timeslice propagators we evaluated the inverse matrix M~^ 
using a single local source on each conhguration, chosen at a site randomly displaced 
from the origin by an even number of lattice spacings in each direction. 

3 The equation of state from RG equations 

Spontaneous breaking of chiral symmetry is signalled by a non-vanishing chiral con¬ 
densate as m —> 0. As already pointed out in the previous section, it is extremely 
difficult to run the HMC for very small values of the bare mass. Therefore we are 
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bound to work with a restricted set of values for m, making an extrapolation to m = 0 
at a fixed value of g rather unreliable. 

In order to hnd numerical evidence of the RG structure described in the intro¬ 
duction, the scaling properties of the theory at the UV hxed point are going to be 
used to hnd an equation of state allowing a global ht of the data. The equation of 
state relates the order parameter for the broken symmetry (in our case the chiral 
condensate) to the value of the external symmetry breaking held (the bare mass m) 
as a function of the external parameters (coupling constant, hnite size of the lattice), 
thus providing a powerful tool for investigating the behaviour of the order parameter 
near the phase transition in the {g,m)-plane. 


3.1 Scaling properties at fixed lattice size 


Let us briehy recall how an equation of state can be obtained from the RG equations 
(RGE) (see eg. [^). Gonsider a theory of a bosonic held (j) regularised with some 
cut-oh A, coupled to an external held H and let M = (0) be the order parameter. 
Then the ehective potential can be expanded in powers of M\ 


r(Af,(,A) = j:— r(")(p, = o,(,A) 


n\ 


(3.1) 


where T^”^ are the IPI Green functions at zero external momenta, t is the reduced 
coupling (see ( 3.10 )) and A the cut-oh. The external held, for given magnetisation, 
can be written as: 


= Ik = ^ (p, = 0,«, A) 


dM ‘—' n\ 

The above equation allows the RGE for H to be written; 

d 


where: 


. d n/ \ 


m = 

git) = 


l + M 


dM 


H{t, m. A) = 0 


A ^ 

-A A log Z 


(3.2) 

(3.3) 

(3.4) 
(3.6) 


and Z is the wave-function renormalisation. 

The solution of the RGE is obtained by using the method of characteristics. In 
a neighbourhood of the hxed point, the solution can be combined with dimensional 
analysis in order to rescale the cut-oh and obtain the general result: 

H(M, t, 1) ~ (3.6) 
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where JF is a universal scaling function. By setting H = 0 in eq. (^), the critical 
behaviour of the order parameter when the external field is switched off is recovered: 


tM ~ const (3.7) 

while, for t = 0, JF(0) is a constant and hence: 

H ^ (3.8) 

showing clearly that (3 and 5 are the usual critical exponents introduced in the context 
of phase transitions. If the critical exponents are related to the existence of a UV 
hxed point, as we are assuming in this section, they must obey the hyper-scaling 
relations, obtained under the assumption that there is a single physical length scale 
near the continuum limit: 

/? = 

5 = 

where rj is the anomalous dimension and v is the critical exponent which characterises 
the divergence of the correlation length as t 0. 

For the Thirring model the reduced coupling t is identified with: 

( = 1/,/ - l/gl (3.10) 


d + 2 


2 + p) 
- 7 ] 


d — 2 + 7] 


(3,9) 


and the symmetry breaking field with the bare mass m, while the order parameter 
is the chiral condensate (xx) defined in the previous section. A Taylor expansion for 
small t reduces Eq. (|T^) to our Erst equation of state: 

m = B{xxY + ^t{xxY~^^^ + O ({t{xx)~^^^f) (3.11) 


If the critical behaviour is described by mean-field theory, the critical exponents 
in Eq. (|3.6| ) are expected to be h = 3 and (3 = 1/2, leading to: 


{xx? 


1 771 At 

B (xx) B 


(3.12) 


In mean-held theory the square of the chiral condensate is a linear function of the 
ratio 7n/{xx)- A plot of {xx? vs. rnl{xx) is called a Fisher plot. In such a plot a 
positive value of the intercept corresponds to a non-vanishing condensate for m = 0, 
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i.e. to chiral symmetry breaking, while the value of the intercept will be exactly zero 
at the critical coupling. 

Eq. ( |3.11| ) permits a five parameter fit of the numerical data for the chiral con¬ 
densate. In order to reduce the number of free parameters in the fit, one can set 
5 — 1/(5 = 1 and keep 5 as the only exponent to be determined from the data. This 
relation originally arose from solution of the SD equations for the gauged Nambu - 
Jona-Lasinio model in four dimensions in ladder approximation here we simply 
use it as a plausible hypothesis to be tested against the data. In what follows we refer 
to the 5 and 4 parameter hts as fits I and II respectively. The outcome of £t I shows 
that the above relation between 5 and (3 is satisfied wihin errors. 

Since both formulae rely on a Taylor expansion around f = 0, they are only 
expected to describe the behaviour of the chiral condensate in a neighbourhood of 
the critical coupling, so that only a reduced set of values of g can be fitted by our 
equation of state. The number of points included in the fit has been chosen in such 
a way to minimise the value. 

In Figs. I and ^ we show the Fisher plots for Nf = 2,4 on the 12^ lattice, together 
with the curves obtained using the parameters of £t II. The existence of a critical 
point separating a chirally symmetric phase from a phase where chiral symmetry is 
spontaneously broken can be easily spotted from the figures according to the above 
criterion. The corresponding lattice data can be fitted to the equation of state, 
yielding stable results and reasonably small The detailed results from the fits are 
reported in Tabs. and |^. Fits I and II agree within errors for both values of on 
the 12^ lattice, justifying the use of the £t with the reduced number of parameters. 
From a practical point of view, such a justification becomes very important since a 
smaller number of variables is more easily fitted from the relatively small datasets that 
one can produce when dealing with dynamical fermions. In the following the results 
quoted will always be the ones from fit II, unless specified. The critical coupling Qc 
does not seem to be unduly affected by finite-size effects, as one can see by comparing 
the results in Tab. B0 The values extracted from the 12^ data for Nj = 2,4 can 
thus be used to check the SD prediction of equation ( |1.9|) : 

£{N, = 2 ) 


9l(N, = 4) 


= 0.291 


Comparing with our result 


a'iiN, - 2) 

9HN, = 4) 


= 0.342 ± 0.016 


( 3 . 13 ) 


( 3 . 14 ) 
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we can conclude that the above prediction is more than three standard deviations 
away from the lattice data. This seems to favour the scenario where a critical number 
of flavors Nf^, exists. For Nf > Nfc the theory would then be chirally symmetric for 
any value of the coupling. This conjecture is also supported by our data for Nf = 6, 
shown in Fig. where it is very difficult to identify a critical point corresponding to a 
symmetry breaking phase transition. Qualitatively as g increases the trajectories seem 
to be accumulating around a line which if continued would intercept the horizontal 
axis. Quantitatively, neither fit I nor fit II provides a satisfactory description of the 
data. It is interesting to note that, in all the cases where the fit converged, the 
critical exponent S is at least two standard deviations away from the mean field value 
5 = 3. Finally, let us remark that finite-size effects influence the value of the critical 
exponent 6: the next sub-section is dedicated to the inclusion of these effects in our 
equation of state, so that the exponents extracted by our fitting procedure are those 
corresponding to the thermodynamic limit. 

3.2 Finite size scaling and the equation of state 

Finite size effects can be used to extract further information about the nature of the 
critical point. In the usual approach to finite size scaling, the RGE are assumed to be 
completely insensitive to finite size effects, since renormalisation deals with the short 
distance properties of the theory; they are therefore unchanged in a finite geometry. 
Nonetheless the flniteness of the volume affects the solutions of RGE, which now 
depend on one more dimensionful parameter. This is equivalent to considering the 
inverse linear size of the lattice, 1/L, as an additional relevant scaling field with 
eigenvalue 1 close to a fixed point at L —>■ oo. By repeating the arguments explained 
above, we obtain an equation expressing the external held as a universal function of 
two rescaled variables: 

MXX)-' ~ (3.15) 

which can again be expanded to yield a new equation of state, 

m = B{xx)^ + -|- higher order terms, (3.16) 

to which the data can be fitted, in principle allowing the extraction of three critical 
exponents (5, (3 and v) and the critical coupling, taking into account finite size effects. 
The parameters obtained from this fit are the thermodynamic (infinite volume) limit 
ones. Hence, assuming the existence of a UV fixed point, the fit to this new equation 
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of state provides a method to determine the thermodynamic limit of the critical 
exponents and coupling without further hypothesis. 

Once again, it is important to stress that Eq. ( |3.16| ) is based on a Taylor expansion 
and is expected to describe the behaviour of the data only in a restricted region around 
the hxed point. As in the previous section, the set of data included in the £t is selected 
by maximising the y^/d.o.f. In order to reduce the number of parameters in the fit, 
(5 —1//3 is set to 1 again and v is assumed to obey the hyperscaling relation ( |3.9| ), as it 
should be in the vicinity of a RG hxed point so that v = (5 + l)/3((5 —1). The number 
of parameters is then reduced from seven to hve. Data from lattice sizes ranging from 
8 to 16 have been included in the £t, whose results are reported in Tab. It is clear 
from the values of the critical exponents that the theory is not mean-held and that 
the data around the critical point seem to match well the expected behaviour of a 
UV hxed point. 

The data for the chiral condensate vs. \l are reported in Fig. and ^ for the 
12^ and 16^ lattices respectively, together with the curves obtained from the ht. The 
dashed line represents the chiral condensate for m = 0 in the thermodynamic limit. 

To conclude, in this section we have shown that the chiral condensate data of 
the Nf = 2 and Nf = 4 models support the hypothesis of a power-law equation of 
state in the critical region. The hts obtained were also consistent with the constraint 
S — 1/13 = 1. In passing it is worth remarking that it is not possible to obtain a 
good ht to our data with a mean-held equation of state with logarithmic corrections, 
as used for QED 4 in [^. The critical exponents extracted for the two models are 
distinct, in neither case resembling those of mean held theory, suggesting that the 
two models have continuum limits described by distinct held theories. Of course, this 
in itself is not direct evidence for an interacting continuum limit at the hxed point; to 
establish that would require the study of lines of constant physics in the space of bare 
parameters g and m, and measurement of n-point functions in the critical region. No 
similar ht was found for Nf = 6, suggesting that in this case there is no transition, 
and that Nfc < 6 for the Thirring model. 


4 Susceptibilities and Spectroscopy 

In this section we review results of susceptibility and spectroscopy measurements, 
principally on the Nf = 2 model studied on a 16^ lattice with m = 0.01 (ie. unless 
stated otherwise), but also showing some results from a smaller lattice, and for Nf = 4, 
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for comparison. First we consider susceptibilities, generically denoted x, which as 
discussed in Sec. 2 are equivalent to bound state propagators evaluated at zero 
three-momentum. They take large values whenever there are long-range correlations 
present. There are two reasons why y is a useful quantity to study. Firstly, it 
can be separated into singlet Xs ^ind non-singlet Xns components; formed from 
disconnected valence fermion lines, is thus a useful measure of the importance of sea 
fermion loops in the vacuum. Since there are no dynamical bosons in the original 
model, any long-range interactions must be mediated by fermion “bubbles”. The 
relative strength of Xs versus Xns is the simplest probe of this effect. Secondly, if we 
assume that fermion - anti-fermion composites resemble canonical bosons, then the 
relation 


X = 


M2 


(4.1) 


holds, where M is the bound state mass and Z is a renormalisation constant associated 
with the projection of the composite operator onto the propagating state. Hence a 
large susceptibility signals a small bound state mass. However, in a non-confining 
theory the assumptions leading to ( |4.1|) may not be justihed, eg. when there is a 
contribution from the fermion - anti-fermion continuum as in the scalar channel of 
the three dimensional GN model P|. In this case extraction of the bound-state mass 
is model dependent, whereas the susceptibility is still unambiguously defined. 

First we consider the scalar channel, plotting in Fig. ^ xis) Xins, and xi versus 
l/g‘^ (see Eq. (|2.25|) ). In the symmetric phase the signal is dominated by xins, which 
peaks at 1 /( 7 ^ ~ 2 — 2.2, indicating a possible light scalar state, before falling away 
in the broken phase. The singlet contribution is almost absent in the symmetric 
phase, but rises across the transition to dominate the signal by l/g"^ = 1.6. The total 
longitudinal susceptibility is approximately constant in the broken phase. 

Next we compare longitudinal and transverse susceptibilities (p.25|,|2.26|) , showing 
data for Nf = 2, m = 0.01 on a 16^ lattice in Fig.|^ and for Nf = A, m = 0.02 on 
a 12^ lattice in Fig. 0. The main features are that xi — Xt deep in the symmetric 
phase, but that Xt rises steeply in the broken phase corresponding to the appearance 
of a massless Goldstone pion. Note that the non-singlet contribution xtns saturates 
the full Ward identity result within errors; this is as expected since singlet - non¬ 
singlet degeneracy could only be lifted by an axial anomaly, which is absent both for 
staggered fermions and for three dimensions. 

There is an interesting relation between the ratio R of longitudinal to transverse 


susceptibilities and the critical scaling behaviour of the model, first noted in [22] and 
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subsequently developed in |^. Using the action (pTT|) and the Ward identity (|2.26|) , 
we have 


2 xi m d{xx) 


Rim,g^) = ^ = 


Xt {xx) dm 

Therefore R can be calculated using the proposed equation of state ( p.ll|) : 

1 


(4.2) 


R{m,g ) = 


+ 


B{XX)' ' 
fjm 


( 4 . 3 ) 


In general R depends on m, but its limiting forms in the chiral limit are easily eval¬ 
uated. In the broken phase, (xx) remains non-vanishing in the chiral limit, whereas 
in the symmetric phase {xx)^ vanishes faster than m; hence 


lim R 

m^O 



9‘^ > gh 
9^ < 9l 


(4.4) 


The near-degeneracy of xi Xt in the symmetric phase from \jg^ = 2.4 onwards 
seen in Fig. 0, and from Xjg^ = 1.2 onwards in Fig. therefore provides additional 
support for our preferred value of^ — l//3 = lin the equation of state (|3.11|) . Exactly 
at the critical coupling, however, R becomes independent of m, and we have 

B\g=gc = (4-5) 

In Tab. ^and Fig. we plot R for various values of m and g near the critical point, in 
particular for four values of m evaluated at our best estimate of the critical coupling 
Xjg^ = 1.92. It is worth noting that the ratio Xis/xins exceeds 40% in this region, 
showing the importance of including the disconnected contributions in this particular 
calculation. However, we also find that R decreases significantly with m, at first sight 
implying that the model is in the symmetric phase at this coupling; the variation of 
R with m is less pronounced aX. 1/g"^ = 1.90, which would thus appear to be closer to 
the critical coupling. This apparent contradiction is resolved if we take the finite-size 
scaling form ( p.l6|) seriously. Using the £t parameters from Tab. || the coefficient 
multiplying the {xxY~^^^ term actually vanishes for 1/g'^ ~ 1.88. We should therefore 
expect to observe a constant R at slightly stronger couplings than the critical one on 
a finite system. Although we have insufficient data to verify this directly, it is clear 
that this kind of analysis has the potential to provide an independent estimate of g^ 
and 5. An optimistic reader would consider the near-agreement found here as a signal 
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of the fact that the RG structure underlying the equation of state approach is indeed 
realised in the Thirring model. 

In Fig. I^we plot the vector susceptibilities Xvs and Xvns ( |2.27|) . The data is quite 
noisy; the trends are that Xvns appears roughly constant in the symmetric phase and 
fall away in the broken phase, and that the ratio Xvs/Xvns lies between 0.2 and 0.3 
throughout the region. 

Next we discuss spectroscopy, which involves the calculation of the timesliced 
correlators C{t) defined in Sec. 2. We follow the usual procedure of extracting the 
mass of the lowest-lying state contributing to the two-point functions from their 
single exponential decay for large time separations. In performing our minimum x^ 
fits, we include correlations between the data at different timeslices. We used the 
bootstrap technique to estimate errors on the fitted parameters. We present results 
from 8^ X 16 and 16^ lattices. 

First consider the fermion, pion and scalar channels. The fermion, since it cannot 
mix with any excited states, shows the clearest signal and can be fitted over the entire 
range by the form 


C,{t) = A , 


( 4 . 6 ) 


where /ur is the physical fermion mass and the minus sign between the forward and 
backward terms is due to our choice of antiperiodic boundary conditions in the time¬ 
like direction. All the quoted results are obtained from fitting in the range t G [2,13] 
with a typical y^/dof ~ 1. Varying the fit range produces negligible variations in the 
fitted parameters and y^/dof. 

Both the scalar and pion channels were fitted by the form 

Ci(t) = A (e-^* -t- . (4.7) 


The pion exhibits single exponential decay from early timeslices and it is possible to 
fit the data in the range t G [3,13] in almost all cases, again with a typical y^/dof ~ 1. 
Similarly to the fermion case, fitting over a smaller time interval within this range 
produces minimal variations. The scalar channel is more problematic. Firstly the 
signals exhibit greater statistical noise. Secondly, as distinct from the pion channel, 
we can only fit the correlators in a restricted time interval, suggesting contamination 
of the signal by higher-mass states. On the 8^ x 16 lattice it is not possible to reliably 
extract the ground-state signal for m = 0.05, while for m = 0.02, fits are possible but 
quite sensitive to the fitted range. In this latter case we therefore choose a conservative 
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fit range t G [4,11]; fitting ontside this range leads to a mnch higher x^/dof. The hts 
on the 16^ data are more stable for all valnes of the fermion bare mass. In the strong 
coupling regime, where the scalar mass gets larger, it proves harder to extract the 
ground-state mass regardless of bare mass or lattice size, resulting in larger errors 
on the htted masses in this phase. In order to improve on this behaviour, it would 
be necessary to move to larger lattices. With regard to this last point, our present 
results should therefore be considered as an attempt to hnd qualitative evidence of 
the phase structure of the model, rather than a precise evaluation of its spectrum. 
The detailed results are presented in Tabs. |^, and Examples of the 

fermion, pion and scalar signals are shown in Figs. |T^ 1^ and ^ respectively. 

A comparison between Tabs. ^ and ^ enables a crude assessment of hnite volume 
effects. In the smaller volume the fermion mass fiR is shifted upwards signihcantly in 
the symmetric phase, while remaining relatively unaffected, within errors, deep in the 
broken phase. The pion, in contrast, is lighter in the smaller volume, the effect again 
being more pronounced in the symmetric phase, but still signihcant in the broken 
phase. The scalar mass seems the least affected by the hnite volume. 

In Fig. 0 we plot the measured masses against l/g"^ for the 16^ system at m = 
0.01. While the pion mass stays constant, the fermion mass is small in the symmetric 
phase, but rises steeply across the transition to become heavier than the pion in the 
broken phase. Just on the symmetric side of the transition we have ~ 2/iR, 
indicating that the pion is a weakly bound state - unfortunately by 1/g^ = 2.4 
this picture breaks down, perhaps due to hnite volume ehects. The scalar mass, 
meanwhile, is just greater but of the same order as the pion mass in the symmetric 
phase, but rises across the transition to become appreciably larger. The scalar and 
pion masses are respectively related to the longitudinal and transverse susceptibilities: 


Xi,i = 


( 4 . 8 ) 


The chiral Ward identity (p^.26|) then yields: 


Mt = 


ixx) 


m. 


( 4 . 9 ) 


As long as chiral symmetry is broken, this equation states that the pion mass squared 
is a linear function of the bare mass and hence vanishes when the breaking is dy¬ 
namical (Goldstone theorem). Fig. shows the behaviour of the pion mass squared 
as m —>■ 0 in the broken phase for Xjg^ = 1.8. The line is a least squares £t; the 
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fitted intercept lies two standard deviations from the origin, which is another signal 
of systematic error. 

It is satisfying that the spectrum reflects the phase transition so sharply, moreover 
in a manner entirely consistent with standard expectations - there is no evidence for 
the unconventional spectrum in the symmetric phase, associated with an infinite order 
phase transition, discussed in []K|. Of course, it would be valuable to refine this initial 
study both by working closer to the chiral limit, and with more lattice spacings in the 
timelike direction to improve the mass estimates - clearly at present we are unable 
to make reliable estimates of binding energies, and hence are unable to make detailed 
comparisons with the predictions of the 1/iVj expansion in the symmetric phase. 

Finally we turn to the vector channels, where we have employed both local ( p.33|) 
and one-link “conserved” ( p.34|) interpolating operators. In either case we found that 
the signal had a strong alternating component, exemplified by the plot shown in 
Fig. Therefore hts of the form ( p.35| ) were used. Superhcially, the main difference 
between the two different vector channels was that the local vector gave numerically 
larger but noisier signals. For both operators we experimented with two parameter 
hts, with B set equal to —A and equal to in (|2.35|) for t G [3,13], then three 
parameter hts in which A and B are allowed to diher, and hnally the most general 
four parameter ht in which A, B, Md and Ma are all independent, the last two forms 
being htted for t G [2,14]. The results are summarised in Tabs. |T^, |T^, and plotted 
in Fig. 

The main trends we found are: 

(i) For all hts the mass of the local operator lies signihcantly above that of the 
conserved operator. This should be no surprise, since as discussed in Sec. 2, and 
displayed in Tab. the two operators have diherent quantum numbers in three 
dimensions. 

(ii) The masses are approximately constant throughout the symmetric phase; across 
the transition the signal becomes smaller and thus noisier, consistent with the falloh 
in Xvns of Fig. |T^. Our mass hts are of limited use for l/g"^ < 1.8. In the conserved 
channel, however, there is evidence for a sharp rise in mass at l/g'^ = 1.9. Hence 
there may be a light vector particle in the symmetric phase, in accordance with the 
l/Nf expansion, but not in the broken phase. 

(in) Although the two parameter hts gave acceptable results, there is tentative evi¬ 
dence in the conserved channel that the four parameter ht is signihcantly better, with 
Ma < Md, particularly from the results with m = 0.02 at l/g"^ = 2.4. An examination 
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of Tab. ^ reveals that the direct channel has the vector spin quantnm nnmber form 
'07/x'0 while the alternating channel has the “axial vector” form Inspection 

of the lattice interaction in the form (|2.12| ) shows that there is a contact interac¬ 
tion between axial vector cnrrents of eqnal strength to that between vector cnrrents; 
therefore it is reasonable that the one-link operator projects onto both states with 
comparable strengths, resnlting in C{t) almost vanishing on even timeslices (Fig. 0). 
Thus we hnd both light vectors and light axial vectors in the symmetric phase. We 
have no explanation for the light axial vector. Following the analysis of P], we can 
examine the pole condition = —M\ in the axial vector channel to leading order in 
l/Nf. 


1 



4m^ — M\ , / 2m -|- Ma \ 

m H-^ In -d. 

AM A \ 2m — Ma / 


0 . 


( 4 . 10 ) 


One readily sees that this equation has no bound state solution of the form Ma = 
2m — £. It would be interesting to see whether this persists at next-to-leading order 
in 1/A^/. 

To summarise, the results in the vector channel present the only unexpected result; 
the light state in the axial vector channel is the first hint of departure from behaviour 
predicted by the l/Nf expansion in the symmetric phase. In future studies at smaller 
bare fermion masses, it would be interesting to monitor the ratio Xvs/Xvns, to see 
whether the channel becomes dominated by fermion bubble chains as in the l/Nf 
expansion. It would also be interesting to repeat this analysis in simulations of non¬ 
compact lattice QEDg to see whether the light axial vector is present or not - in this 
case since the held is now a dynamical variable, it is possible that short-wavelength 
huctuations will be suppressed by a factor 1/fc^, preventing the strong coupling to 
the axial vector channel occurring in (|2.12|) . 


5 Conclusions 

In conclusion we briehy summarise our main hndings, and suggest further directions 
to explore. For Nf = 2 and Nf = 4 the Thirring model exhibits a phase transition 
at hnite critical coupling g/ to a phase in which chiral symmetry is spontaneously 
broken. This contradicts the naive l/Nf expansion, but supports the predictions of 
non-perturbative methods such as the Schwinger-Dyson approach. In the vicinity of 
the critical point the model’s equation of state seems well htted by a renormalisation- 
group inspired form corresponding to power-law scaling. The picture remains con- 
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sistent for Nf = 2 even once a finite volnme scaling analysis, nsing data from larger 
lattices and smaller bare masses, is employed. The critical exponents we find differ 
significantly between Nf = 2 and Nf = 4, snggesting that the two models lie in 
different nniversality classes. Althongh it is tempting to infer the existence of non- 
perturbative continnnm limits at these points, it would be nice to have direct evidence 
by measurement of a renormalised coupling constant via a 3- or 4-point correlation 
function. 

For Nf = 6 the picture is not so clear, and we did not hnd a hxed-point £t for the 
equation of state. The Fisher plot of Fig. strongly suggests that in the chiral limit 
the model lies in the symmetric phase for all values of the coupling, and hence that 
Nfc < 6. However, bearing in mind the practical difficulties in identifying the strong 
coupling limit discussed in Sec. 2, perhaps a more conservative conclusion is that the 
Nf = 6 model seems qualitatively very different from Nf = 2,4 It might be 

very difficult ever to exclude absolutely a chiral condensate exponentially suppressed 
in Nf for large Nf] however our data disagrees with the only theoretical prediction 
( P-.QI) to put forward this scenario [^. We intend to extend the Nf = 6 studies to 
larger lattices and smaller m to pursue this issue further. 

Although we made no attempt to test the scaling form ( |1.3| ), corresponding to 
an inhnite order phase transition, our studies of the susceptibilities and spectrum 
prejudice us against this scenario: the relative ordering of the physical fermion, scalar 
and pseudoscalar masses varies continuously across the transition; there appear to be 
light bound states in the symmetric phase, comparable with the scale set by the 
physical fermion mass, in contradiction to the behaviour expected near a conformal 
transition discussed in Refs. [|^, 0. It would be interesting to extend this study to 
Nf = 4, where Qc lies closer to the strong coupling limit, to see if the spectrum shows 
any qualitative difference. Work in this direction is in progress. 

Finally, although our admittedly exploratory spectrum analysis has yielded im¬ 
portant information, it is clear that larger lattices in the timelike direction will be 
needed before accurate measurements of masses and binding energies, or indeed de¬ 
partures from the exponential form of the decay corresponding to a spectral function 
more complicated than simple isolated poles, can be accomplished. Accurate mass 
measurements would be desirable in order to compare with predictions from the 1/Nf 
expansion, to check whether the latter approach has any validity near the transition 
even in the symmetric phase. There is a first hint of its breaking down in the appar¬ 
ent existence of a light axial vector state here; however it is clear much more work is 





needed both numerically and analytically. 
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Table 1: List of results for the chiral condensate for Nf = 2 


L 

m 


(xx) 

A(xx) 

8 

0.01 

1.8 

0.08610 

0.0030 

8 

0.01 

2.0 

0.06130 

0.0020 

8 

0.01 

2.2 

0.05170 

0.0020 

8 

0.01 

2.4 

0.00425 

0.0010 

8 

0.02 

0.5 

0.25097 

0.0084 

8 

0.02 

1.0 

0.28171 

0.0110 

8 

0.02 

1.5 

0.20172 

0.0068 

8 

0.05 

1.6 

0.26956 

0.0080 

8 

0.05 

1.8 

0.22639 

0.0060 

8 

0.05 

2.0 

0.20417 

0.0040 

8 

0.05 

2.2 

0.18051 

0.0043 

8 

0.05 

2.4 

0.15785 

0.0038 

8 

0.05 

2.6 

0.14403 

0.0035 

8 

0.05 

2.8 

0.13393 

0.0028 

12 

0.01 

1.6 

0.17485 

0.0057 

12 

0.02 

1.6 

0.22341 

0.0055 

12 

0.02 

1.8 

0.19083 

0.0047 

12 

0.02 

2.0 

0.14680 


12 

0.02 

2.2 

0.11172 

0.0026 

12 

0.02 

2.4 

0.08992 


12 

0.02 

2.6 

0.07737 

0.0015 

12 

0.02 

2.8 

0.06493 


12 

0.02 

3.0 

0.05715 


12 

0.02 

3.2 

0.05256 


12 

0.02 

3.4 

0.04784 


12 

0.03 

1.6 

0.23970 

0.0035 

12 

0.03 

1.8 

0.20873 

0.0036 

12 

0.03 

2.0 

0.17070 

0.0023 

12 

0.03 

2.2 

0.14342 

0.0025 

12 

0.03 

2.4 

0.12171 

0.0017 

12 

0.03 

2.6 

0.10508 

0.0013 

12 

0.03 

2.8 

0.09153 

0.0011 

12 

0.03 

3.0 

0.08303 


12 

0.03 

3.2 

0.07554 


12 

0.03 

3.4 

0.07107 
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Table 2: List of results for Nf = 2 (continued) 


L 

m 

i/g^ 

(xx) 

A(xx) 

12 

0.04 

1.6 

0.25599 

0.0031 

12 

0.04 

1.8 

0.22321 

0.0032 

12 

0.04 

2.0 

0.19460 

0.0030 

12 

0.04 

2.2 

0.17244 

0.0024 

12 

0.04 

2.4 

0.14307 

0.0019 

12 

0.04 

2.6 

0.12966 

0.0018 

12 

0.04 

2.8 

0.11604 

0.0015 

12 

0.04 

3.0 

0.10851 

0.0012 

12 

0.04 

3.2 

0.09707 

0.0011 

12 

0.04 

3.4 

0.09282 

0.0008 

12 

0.05 

1.6 

0.27098 

0.0030 

12 

0.05 

1.8 

0.23624 

0.0028 

12 

0.05 

2.0 

0.20952 

0.0025 

12 

0.05 

2.2 

0.18450 

0.0020 

12 

0.05 

2.4 

0.16750 

0.0020 

12 

0.05 

2.6 

0.15171 

0.0020 

12 

0.05 

2.8 

0.13861 

0.0014 

12 

0.05 

3.0 

0.12572 

0.0012 

12 

0.05 

3.2 

0.11861 

0.0011 

12 

0.05 

3.4 

0.11118 

0.0010 

16 

0.01 

1.6 

0.20250 

0.0028 

16 

0.01 

1.8 

0.14190 

0.0021 

16 

0.01 

1.9 

0.12780 

0.0032 

16 

0.01 

1.92 

0.12419 

0.0018 

16 

0.01 

2.0 

0.10250 

0.0021 

16 

0.01 

2.1 

0.08614 

0.0020 

16 

0.01 

2.2 

0.07177 

0.0012 

16 

0.01 

2.4 

0.05287 

0.0010 

16 

0.02 

1.6 

0.22040 

0.0030 

16 

0.02 

1.92 

0.15733 

0.0013 

16 

0.02 

2.0 

0.14580 

0.0012 

16 

0.02 

2.4 

0.09430 

0.0006 

16 

0.03 

1.8 

0.20440 

0.0010 

16 

0.03 

1.9 

0.18624 

0.0009 

16 

0.03 

1.92 

0.18547 

0.0010 

16 

0.03 

2.0 

0.17174 

0.0009 

16 

0.03 

2.1 

0.15632 

0.0009 

16 

0.04 

1.92 

0.20576 

0.0008 
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Table 3: 

List of results for Nf = 4 


L 

m 

i/g^ 

ixx) 

Mxx) 

12 

0.02 

0.6 

0.21771 

0.0037 

12 

0.02 

0.7 

0.17435 

0.0032 

12 

0.02 

0.8 

0.13163 


12 

0.02 

0.9 

0.09994 

0.0024 

12 

0.02 

1.0 

0.07917 

0.0018 

12 

0.02 

1.1 

0.06700 

0.0011 

12 

0.02 

1.2 

0.06030 


12 

0.02 

1.3 

0.05453 


12 

0.02 

1.4 

0.04882 


12 

0.02 

1.5 

0.04636 


12 

0.03 

0.6 

0.23476 

0.0033 

12 

0.03 

0.7 

0.20343 

0.0031 

12 

0.03 

0.8 

0.16221 

0.0028 

12 

0.03 

0.9 

0.13419 


12 

0.03 

1.0 

0.11065 

0.0017 

12 

0.03 

1.1 

0.09750 


12 

0.03 

1.2 

0.08977 

0.0016 

12 

0.03 

1.3 

0.07804 


12 

0.03 

1.4 

0.07295 


12 

0.03 

1.5 

0.06793 


12 

0.04 

0.6 

0.24741 

0.0036 

12 

0.04 

0.7 

0.21080 

0.0028 

12 

0.04 

0.8 

0.18680 

0.0024 

12 

0.04 

0.9 

0.16079 


12 

0.04 

1.0 

0.13808 


12 

0.04 

1.1 

0.12140 

0.0018 

12 

0.04 

1.2 

0.10850 

0.0013 

12 

0.04 

1.3 

0.10094 


12 

0.04 

1.4 

0.09362 


12 

0.04 

1.5 

0.08734 


12 

0.05 

0.6 

0.25489 

0.0026 

12 

0.05 

0.7 

0.23618 

0.0025 

12 

0.05 

0.8 

0.20854 

0.0022 

12 

0.05 

0.9 

0.18000 

0.0021 

12 

0.05 

1.0 

0.160 

0.0016 

12 

0.05 

1.1 

0.14423 

0.0015 

12 

0.05 

1.2 

0.12921 

0.0013 

12 

0.05 

1.3 

0.12051 

0.0012 

12 

0.05 

1.4 

0.11202 

0.0011 

12 

0.05 

1.5 

0.10653 
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Table 4: List of results for iV/ = 6 


L 

m 

i/g^ 

(xx) 

A(xx) 

12 

0.02 

0.25 

0.14342 

0.0034 

12 

0.02 

0.3 

0.14687 

0.0035 

12 

0.02 

0.35 

0.13924 

0.0040 

12 

0.02 

0.4 

0.11832 

0.0030 

12 

0.02 

0.45 

0.09867 

0.0040 

12 

0.02 

0.5 

0.07608 

0.0020 

12 

0.03 

0.25 

0.16721 

0.0030 

12 

0.03 

0.3 

0.17836 

0.0040 

12 

0.03 

0.35 

0.17099 

0.0032 

12 

0.03 

0.4 

0.15218 

0.0028 

12 

0.03 

0.45 

0.12766 

0.0030 

12 

0.03 

0.5 

0.10716 

0.0020 

12 

0.04 

0.3 

0.18403 

0.0029 

12 

0.04 

0.35 

0.18659 

0.0031 

12 

0.04 

0.4 

0.17549 

0.0026 

12 

0.04 

0.45 

0.15650 

0.0022 

12 

0.04 

0.5 

0.13959 

0.0038 

12 

0.05 

0.25 

0.19691 

0.0028 

12 

0.05 

0.3 

0.20373 

0.0029 

12 

0.05 

0.35 

0.20636 

0.0028 

12 

0.05 

0.4 

0.19784 

0.0030 

12 

0.05 

0.45 

0.17549 

0.0017 

12 

0.05 

0.5 

0.15903 

0.0020 


Table 5: Spin/flavor 

assignments of the fermion bilinears used in spectroscopy 


direct 

alternating 

pion 

75 ® 1 


scalar 

1 (g) 1 

® r*T* 

local vector 

7/.73 ® T-gV; 

757/.73 ® r^r* 

conserved vector 

7/. ® 1 

757m ® ^3 


35 










Table 6: Results from fits on the 12^ lattice. 



Parameter 

Fit I 

Fit II 

Nf = 2 

ygi 

2.03(9) 

1.94(4) 


5 

2.32(23) 

2.68(16) 


P 

0.71(9) 

— 


A 

0.32(5) 

0.37(1) 


B 

1.91(43) 

2.86(35) 


xVd-o.f 

2.4 

2.1 

Nf = 4 

i/gl 

0.63(1) 

0.66(1) 


6 

3.67(28) 

3.43(19) 


P 

0.38(4) 

— 


A 

0.78(5) 

0.73(2) 


B 

7.9(2.8) 

6.4(1.5) 


xVd.o.f 

3.1 

2.0 


Table 7: Results from fits on the 16^ lattice. 


Parameter 

Fit I Fit II 

Nf = 2 

yg'i 

1.93(4) 


5 

2.55(15) 


P 

— 


A 

0.38(1) 


B 

2.29(35) 


xVd-o.f 

2.3 


Table 8: Results from fit including finite 

size scaling. 


Parameter 

Fit III 

Nf = 2 

yg'i 

1.92(2) 


6 

2.75(9) 


Pn 

0.57(2) 


vh 

0.60(2) 



0.71(4) 


A 

0.334(7) 


B 

2.7(3) 


C 

2.1(7) 


xVd.o.f 

1.76 


“evaluated from 5 — \/(5 = 1 constraint 
^evaluated from hyperscaling relation 
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Table 9: Longitudinal and transverse susceptibilities in the vicinity of the critical 
coupling. 



m 

Xls 

Xlns 

Xi 

Xt R = xi/xt 

1.9 

0.01 

1.51(15) 

3.40(40) 

4.90(43) 



0.03 

0.57(8) 

1.68(15) 

2.24(17) 


1.92 

0.01 

1.65(19) 

3.88(28) 

5.53(34) 

12.17(25) 0.454(29) 


0.02 

0.72(8) 

2.41(15) 

3.13(17) 

7.83(7) 0.400(22) 


0.03 

0.67(9) 

1.70(10) 

2.37(13) 

6.17(5) 0.384(21) 


0.04 

0.34(5) 

1.35(8) 

1.69(9) 

5.15(3) 0.329(18) 

2.0 

0.01 

1.13(19) 

4.25(44) 

5.37(48) 

10.12(21) 0.531(49) 


0.03 

0.62(8) 

1.93(9) 

2.54(12) 

5.71(5) 0.445(21) 


Table 10 

Masses for the fermion, pion and scalar, 16^ lattice, m = 0.01 

1/17^ 


hi? 



Ms 

1.6 


0.36 (3) 


0.20 (1) 

— 

1.8 


0.23 (2) 


0.23 (2) 

0.47 (4) 

1.9 


0.16 (1) 


0.21 (1) 

0.33 (2) 

2.0 


0.10 (1) 


0.24 (1) 

0.29 (1) 

2.1 


0.10 (1) 


0.22 (1) 

0.27 (1) 

2.2 


0.08 (1) 


0.21 (1) 

0.26 (1) 

2.4 


0.06 (1) 


0.21 (1) 

0.23 (1) 


Table 11 

Masses for the fermion, pion and scalar, 16^ lattice, m = 0.02 

1/77^ 


fJ'R 



Ms 

1.6 


0.330 (40) 


0.280 (4) 

0.52 (15) 

1.8 


0.360 (40) 


0.290 (6) 

0.39 (8) 

1.92 


0.200 (10) 


0.280 (4) 

0.43 (6) 

2.0 


0.190 (10) 


0.300 (5) 

0.35 (5) 

2.4 


0.109 (3) 


0.290 (4) 

0.27 (1) 


Table 12 

Masses for the fermion, pion and scalar, 16^ lattice, m = 0.03 

1/77^ 


I^R 



Ms 

1.8 


0.290 (20) 


0.340 (5) 

0.64 (9) 

1.9 


0.240 (20) 


0.338 (4) 

0.60 (4) 

1.92 


0.250 (20) 


0.338 (4) 

0.64 (6) 

2.0 


0.270 (20) 


0.344 (4) 

0.54 (4) 

2.1 


0.210 (10) 


0.341 (4) 

0.50 (3) 
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Table 13: Masses for the fermion, pion and scalar, 8^ x 16 lattice, m = 0.02 


l/f 

hr? 


Ms 

1.6 

**** 


**** 

1.7 

0.340 (30) 

0.256 (4) 

0.42 (5) 

1.9 

0.290 (20) 

0.240 (5) 

0.44 (7) 

2.0 

0.280 (10) 

0.231 (6) 

0.30 (7) 

2.1 

0.200 (10) 

0.231 (6) 

0.26 (3) 

2.2 

0.240 (10) 

0.225 (5) 

0.29 (2) 

2.3 

0.185 (5) 

0.225 (6) 

0.23 (2) 

2.4 

0.170 (9) 

0.206 (5) 

0.26 (2) 


Table 14: Masses for the fermion, pion and scalar, 8^ x 16 lattice, m = 0.05 


i/g'^ 

l^R 


Ms 

1.6 

0.480 (40) 

0.398 (4) 

— 

1.7 

0.470 (30) 

0.390 (4) 

— 

1.8 

— 

0.390 (4) 

— 

1.9 

0.350 (10) 

0.386 (3) 

— 

2.0 

0.350 (20) 

0.378 (4) 

— 

2.1 

0.320 (10) 

0.368 (4) 

— 

2.2 

0.280 (10) 

0.362 (4) 

— 

2.3 

— 

0.359 (3) 

— 

2.4 

0.255 (9) 

0.349 (3) 

— 

2.5 

0.246 (7) 

0.340 (3) 

— 
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Table 15: Masses obtained in the local vector channel from a 16^ lattice 


htted parameters 

m 

i/g'^ 

Md 

Ma 

xVdof 

2 

0.01 

3.0 

0.307(18) 


7.1 



2.4 

0.332(63) 


1.4 



2.2 

0.357(79) 


4.2 



2.1 

0.379(109) 


1.2 



2.0 

0.360(61) 


1.7 



1.9 

0.389(88) 


0.4 



1.8 

0.807(665) 


2.3 


0.02 

2.4 

0.365(34) 


4.2 


0.03 

2.1 

0.952(183) 


3.7 



2.0 

0.870(171) 


2.6 



1.9 

0.438(88) 


2.5 



1.8 

0.928(333) 


1.5 

3 

0.01 

3.0 

0.310(19) 


3.4 



2.4 

0.350(57) 


0.8 



2.2 

0.388(67) 


1.8 



2.1 

0.479(132) 


1.5 



2.0 

0.358(59) 


1.0 



1.9 

0.397(87) 


0.3 



1.8 

0.764(495) 


2.4 

4 

0.01 

3.0 

0.263(18) 

0.256(19) 

1.2 



2.4 

0.309(57) 

0.294(62) 

0.7 



2.2 

0.336(66) 

0.284(100) 

1.9 



2.1 

0.448(106) 

0.312(103) 

1.1 



2.0 

0.358(60) 

0.345(79) 

1.3 



1.9 

0.395(100) 

0.337(173) 

0.5 



1.8 

0.701(457) 

1.368(2.7) 

2.5 
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Table 16: Masses obtained in the conserved vector channel from a 16^ lattice 


htted parameters 

m 

i/g^ 

Md 

Ma 

xVdof 

2 

0.01 

3.0 

0.237(18) 


3.0 



2.4 

0.259(30) 


2.1 



2.2 

0.255(29) 


1.3 



2.1 

0.271(45) 


0.7 



2.0 

0.244(53) 


1.5 



1.9 

0.386(80) 


0.5 



1.8 

0.409(224) 


0.6 


0.02 

2.4 

0.306(18) 


2.2 


0.03 

2.1 

0.489(45) 


1.3 



2.0 

0.656(75) 


1.9 



1.9 

0.566(114) 


2.8 



1.8 

0.814(218) 


0.9 

3 

0.01 

3.0 

0.239(15) 


2.8 



2.4 

0.266(30) 


1.9 



2.2 

0.257(29) 


1.5 



2.1 

0.312(42) 


0.5 



2.0 

0.258(56) 


1.8 



1.9 

0.456(89) 


0.9 



1.8 

0.974(879) 


0.9 

4 

0.01 

3.0 

0.253(18) 

0.230(17) 

2.9 



2.4 

0.293(34) 

0.229(36) 

1.8 



2.2 

0.325(37) 

0.162(46) 

0.5 



2.1 

0.320(43) 

0.259(73) 

0.5 



2.0 

0.302(63) 

0.223(67) 

1.9 



1.9 

0.605(204) 

0.190(203) 

0.5 


0.02 

2.4 

0.366(20) 

0.220(23) 

1.2 
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Figure 1: /9-function for d = 2 and 2 < d < 4, showing for the latter case a non-trivial 
UV hxed point and phase structure. 



Figure 2: Diagrams contributing to vacuum polarisation in lattice QED. 






Figure 5: Fisher plot for Nf = 4 on 12^ lattice. Values of 1/g^ range from 0.5 (Q) 
steps of 0.1 to 1.4 (*). 



Figure 6; Fisher plot for Nf = 6 on 12^ lattice. 





0.40 



Figure 7: Chiral condensate vs. l/g"^ on 12^ lattice. The solid lines are the hts from 
(p.l6|), the dashed line the extrapolation to chiral and thermodynamic limits. 



Figure 8: Chiral condensate vs. 1/g'^ on 16^ lattice. 
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Figure 11: Xt (squares), Xtns (triangles) and xins (circles) vs. 1/g^ for Nj = 4. 



Figure 12: Susceptibility ratio R vs. m for 1/g^ = 1.9 (circles), 1.92 (squares) and 
2.0 (diamonds). The horizontal band marks the value of l/h from the fit (p.l6|) . 
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Figure 13: Vector susceptibilities Xvs (circles) and Xvns (diamonds) vs. 1/g 



Figure 14: Fermion correlator for 1/g'^ = 2.2. 















